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Abstract. We show that quatcrnionic Gaussian random variables satisfy a 
generalization of the Wick formula for computing the expected value of prod- 
ucts in terms of a family of graphical enumeration problems. When applied to 
the quatcrnionic Wigncr and Wishart families of random matrices the result 
gives the duality between moments of these families and the corresponding 
real Wigncr and Wishart families. 



1. Introduction 

The duality between symplectic and orthogonal groups has a long standing his- 
tory, and has been noted in physics literature in various settings, see e.g. [151 EQI 
\T7\ . Informally, the duality asserts that averages such as moments or parti- 
tion functions for the symplectic case of "dimension" TV, can be derived from the 
respective formulas for the orthogonal case of dimension N by inserting —N into 
these expressions and by simple scaling. The detailed study of the moments of 
one-matrix Wishart ensembles, with duality explicitly noted, appears in |10j . see 
[ini Corollary 4.2]. The duality for one matrix Gaussian Symplectic Ensemble was 
noted by Mulase and Waldron [21] who introduced Mobius graphs to write the ex- 
pansion for traces of powers of GOE/GUE/GSE expansions in a unified way. The 
duality appears also in [17l Theorem 6] as a by-product of differential equations 
for the generating functions of moments. Ref. [3 [71 [8] analyze the related "genus 
series" over locally orientable surfaces. 

The purpose of this paper is to prove that the duality between moments of 
the Gaussian Symplectic Ensemble and the Gaussian Orthogonal Ensemble, and 
between real Wishart and quatcrnionic Wishart ensembles extends to several inde- 
pendent matrices. Our technique consists of elementary combinatorics; our proofs 
differ from |21| in the one matrix case, and provide a more geometric interpre- 
tation for the duality; in the one-matrix Wishart case, our proof completes the 
combinatorial approach initiated in |10[ Section 6] . The technique limits the scope 
of our results to moments, but the relations between moments suggest similar re- 
lations between other analytic objects, such as partition functions, see [21], [16]. 
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The asymptotic expansion of the partition function and analytic description of the 
coefficients of this expansion for /3 = 2 case appear in [H [9l [19] . 

The paper is organized as follows. In Section [2] we review basic properties 
of quatcrnionic Gaussian random variables. In Section [3] wc introduce Mobius 
graphs; Theorems 13.11 and 13.21 give formulae for the expected values of products 
of quatcrnionic Gaussian random variables in terms of the Euler characteristics of 
sub- families of Mobius graphs or of bipartite Mobius graphs. In Section [4] we apply 
the formulae to the quatcrnionic Wigner and Wishart families. 

In this paper, we do not address the question of whether the duality can be 
extended to more general functions, or to more general /3-Hermite and /3-Laguerre 
ensembles introduced in [3]. 



2. Moments of quaternion-valued Gaussian random variables 

2.1. Quaternion Gaussian law. Recall that a quaternion q G M can be repre- 
sented as q = + ixi + jx2 + kx^ with = p = = ijk = — 1 and with real 
coefficients xq, . . . ,2^3. The conjugate quaternion is 9 = xo — ixi — jx2 — kx^, so 
9p 99 ^ 0. Quaternions with xi = X2 = X3 ^ are usually identified with real 
numbers; the real part of a quaternion is ^{q) ~ {q + 9)/2. 

It is well known that quaternions can be identified with the set of certain 2x2 
complex matrices: 



(2.1) M3 xo+ ixi + jx2 + kxz 



.To -I- IXi X2 + 1x3 

-X2 + ixz xq ~ ixi 



e A^2x2(C), 



where on the right hand side i is the usual imaginary unit of C. Note that since 
5R((?) is twice the trace of the matrix representation in (|2.ip . this implies the cyclic 
property 

(2.2) 3?(q'i(72 . . . 9n) = 5?(g2g3 ■ • ■ Qriqi)- 

The (standard) quaternion Gaussian random variable is an H-valued random 
variable which can be represented as 

(2.3) Z = ^„+iCi+j^2 + k^3 

with independent real normal iV(0, 1) random variables Co7CiiC2,?3- Due to sym- 
metry of the centered normal laws on R, the law of (Z, Z) is the same as the law 
of {Z,Z). A calculation shows that if Z is quaternion Gaussian then for fixed 

qi,q2 e H, 

E{ZqiZq2) = E{Z^)qiq2, E{Zq{Zq2) = E{ZZ)^{qi)q2 . 
For future reference, we insert explicitly the moments: 

(2.4) nZqiZq2) = -2gi(72, 

(2.5) E{ZqiZq2) = 2{qi+qi)q2. 

By linearity, these formulas imply 

(2.6) E(5R(Zgi)»(Zg2)) = 5R(<?ig2), 

(2.7) E{^{Zq,MZq2)) = miq2)- 
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2.2. Moments. The following is known as Wick's theorem [26] . 

Theorem A (Isscrlis US])- If {Xi, . . . , X2n) is a M."^"" -valued Gaussian random 
vector with mean zero, then 

(2.8) EiXiX2...X2n) n Ei^J^k), 

V {j,k}&V 

where the sum is taken over all pair partitions V of {1, 2, . . . , 2n}, i.e., partitions 
into two-element sets, so each V has the form 

V = {{jl, fcl}, {j2, fc2}, • ■ • , {jri, fcn}} • 

Theorem [X] is a consequence of the moments-cumulants relation [18] ; the con- 
nection is best visible in the partition formulation of [23. For another proof, see 
[H page 12]. 

Our first goal is to extend this formula to certain quaternion Gaussian ran- 
dom variables. The general multivariate quaternion Gaussian law is discussed in 
[24| . Here we will only consider a special setting of sequences that are drawn with 
repetition from a sequence of independent standard Gaussian quaternion random 
variables. In section [4] we apply this result to a multi-matrix version of the duality 
between GOE and GSE ensembles of random matrices. 

In view of the Wick formula (j2.8p for real-valued jointly Gaussian random vari- 
ables, formulas (|2.4[) and (|2.5p allow us to compute moments of certain products 
of quaternion Gaussian random variables. Suppose the n-tuple {Xi,X2, . . . ,^n) 
consists of random variables taken, possibly with repetition, from the set 

{Zi, Zi, Z2, Z2, • . ■ }, 

where Zi, Z2, . . . arc independent quaternion Gaussian random variables. Consider 
an auxiliary family of independent pairs { (^^'"\ ^^^'"^ ) : r = 1,2, . . .} which have 
the same laws as (Xj, X^), 1 < j, k < n and are independent for different r. Then 
the Wick formula for real-valued Gaussian variables implies ¥.{XiX2 . . . Xn) = 
for odd n, and 

(2.9) E{XiX2 ...Xn) = ^E(y/^(^^V2^^^^^) . . .Fi-^^")^), 

/ 

where the sum is over the pair partitions V that appear under the sum in Theorem 
\M each represented by the level sets of a two-to-one valued function / : {1, . . . , n} — s- 
{1, . . . , m} for n = 2m. (Thus the sum is over classes of equivalence of /, each of 
ml representatives contributing the same value.) 

For example, if Z is quaternion Gaussian then applying (j2.9|) with /i that is 
constant, say 1, on {1,2}, /2 that is constant on {1,3}, and f^ that is constant on 
{1, 4} we get 

E{Z'^) = E{Z^) {E{Z^) + E{ZZ) + E{Z^)) = 0. 

Formulas (|2.4p and (|2.5p then show that the Wick reduction step takes the fol- 
lowing form. 

(2.10) E{XiX2 . ..X„) = ^E(XiX,)E(C/,X,+i . . .X„), 

J=2 
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where 

Uj = <Xj^i...X2 iiXj=Xi 
[ otherwise . 

This imphes that one can do inductively the calculations, but due to noncommu- 
tativity arriving at explicit answers may still require significant work. 

Formula (j2.10p implies that £(^1X2 . . . X„) is real, so on the left hand side of 
()2.10p we can write E(5ft(XiX2 . . . Xn)); this form of the formula will be associated 
with one- vertex Mobius graphs. 

Furthermore, we have a Wick reduction which will correspond to the multiple 
vertex Mobius graphs: 

(2.11) E(3?(Xi)5R(X2X3...X„)) 

n 

= ^E(5R(Xi)5R(X,))E(5R(X2 . . . Xj-i^j+i . • . X„)). 

(This is just a consequence of Theorem IX| . 

In the next section we will show that formulae (|2.10|1 and (j2.11|l give a method 
of computing the expected values of quaternionic Gaussian random variables by 
the enumeration of Mobius graphs partitioned by their Euler characteristic. This 
is analogous to similar results for complex Gaussian random variables and ribbon 
graphs, and for real Gaussian random variables and Mobius graphs. In Section [4] 
we will show that this result implies the duality of the GOE and GSE ensembles of 
Wigner random matrices and the duality of real and quaternionic Wishart random 
matrices. 

3. Mobius graphs and quaternionic Gaussian moments 

In this section we introduce Mobius graphs and then give formulae for the ex- 
pected values of products of quaternionic Gaussian random variables in terms of 
the Euler characteristics of sub- families of Mobius graphs. This is an analogue of 
the method of t'Hooft [U [5l [H (TTJ [131 US]. Mobius graphs have also been used 
to give combinatoric interpretations of the expected values of traces of Gaussian 
orthogonal ensemble of random matrices and of Gaussian symplectic ensembles, see 
the articles [5] and [U]. The connection between Mobius graphs and quaternionic 
Gaussian random variables is at the center of the work of Mulase and Waldron [5T] . 

3.1. Mobius graphs. Mobius graphs are ribbon graphs where the edges (ribbons) 
are allowed to twist, that is they either preserve or reverse the local orientations of 
the vertices. As the convention is that the ribbons in ribbon graphs are not twisted 
we follow [21] and call the unoriented variety Mobius graphs. The vertices of a 
Mobius graph are represented as disks together with a local orientation; the edges 
are represented as ribbons, which preserve or reverse the local orientations of the 
vertices connected by that edge. Next we identify the collection of disjoint cycles 
of the sides of the ribbons found by following the sides of the ribbon and obeying 
the local orientations at each vertex. We then attach disks to each of these cycles 
by gluing the boundaries of the disks to the sides of the ribbons in the cycle. These 
disks are called the faces of the Mobius graph, and the resulting surface we find 
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is the surface of maximal Euler characteristic on which the Mobius graph may be 
drawn so that edges do not cross. 

Denote by f (F), e(r), and /(F) the number of vertices, edges, and faces of F. 
We say that the Euler characteristic of F is 

x(F) = t;(F)-e(r) + /(r), 

for connected F, this is also the maximal Euler characteristic of a connected surface 
into which F is embedded. For example, in Fig. [l] the Euler characteristics are 
Xi = l — 1 + 2 = 2 and X2 = l — 1 + 1 = 1- The two graphs may be embedded into 
the sphere or projective sphere respectively. 

If F decomposes into connected components Fi, F2, then x(F) = x(ri) + x(F2)- 
Throughout the paper our Mobius graphs will have the following labels attached 
to them: the vertices are labeled to make them distinct, in addition the edges 
emanating from each vertex are also labeled so that rotating any vertex produces 
a distinct graph. These labels may be removed if one wishes by rescaling all of our 
quantities by the number of automorphisms the unlabeled graph would have. 

3.2. Quaternion version of Wick's theorem. Suppose the 2n-tuple 

{Xi, X2, . . . , X2n) 

consists of random variables taken, possibly with repetition, from the set 
{Zi, Zi, Z2, Z2, ■ ■ ■} where Zi, Z2, ■ ■ ■ are independent quaternionic Gaussian. Fix 
a sequence ji , ^2 , • ■ • , jm of natural numbers such that ji + ■ ■ ■ + jm = 2n. 

Consider the family M = Mji^...j^{Xi,X2, ■ ■ ■ ,X2n), possibly empty, of 
Mobius graphs with m vertices of degrees ji,j2,---,jm with edges labeled by 
Xi, X2, . . . , X2„, whose regular edges correspond to pairs Xi = Xj and flipped 
edges correspond to pairs Xi = Xj. No edges of F G can join random variables 
Xi,Xj that are independent. 

Theorem 3.1. Let {Xi,X2, ■ ■ ■ , X2n} be chosen, possibly with repetition, from the 
set {Zi, Zi, Z2, Z2, ■ ■ ■} where Zj are independent quaternionic Gaussian random 
variables, then 

(3.1) E(5R(XiX2 . ..X,,MX^^^, . . X . . . 

X 3'?(-''^ji+j2 + ...+j„_i+l • • • X2n)) 

(The right hand side is interpreted as when M ~ ib.) 

Remark 3.1. We would like to emphasize that in computing the Euler characteristic 
one must first break the graph into connected components. For example, if ji = 
• • • = j„i = 1 so that m = 2?! is even, and Xi,X2, ■ ■ ■ ,X2n are n independent 
pairs, as the real parts are commutative we may assume that X2k = ^2fc-i, and 
the moment is 

1 = E(3fJ(Xi)5R(X2) . . . n{X2n)) = 4-"(~2)>^(r)^ 

We see that graphically F is a collection of 2n degree one vertices connected together 
forming n dipoles (an edge with a vertex on either end) . Hence there are n connected 
components each of Euler characteristic 2, therefore the total Euler characteristic 
is X = 2n giving 

4-"{-2)^ = 4""4" = 1. 



E(-2) 



x(r) 
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Proof of Theorem\3Jl In view of l[2?9l) and fTTOl) . it suffices to sliow that if 
Xi, . . . , X2n consists of n independent pairs, and each pair is either of the form 
{X,X) or {X,X), then 

(3.2) E(5R(XiX2 . ..X,,MX^^^, . . X . . . 

X . ..X,J) = 4"-™(-2)>^(r) 

where F is the Mobius graph that describes the pairings of the sequence. 
First we check the two Mobius graphs for 71 = 1, m = 1: 

E{n{XX)) = (-2)2 , and E(3?(XX)) = (-2)^ . 

One checks that these correspond to the Mobius graphs in Figure [1] which gives a 
sphere (x = 2) and projective sphere (x ~ ^) respectively. 




XX XX 

Figure 1 . The two possible Mobius graphs with a single degree 2 
vertex. The left hand one is a ribbon which is untwisted and the 
graph embeds into a copy of the Ricmann sphere, while the right 
hand one is a ribbon which is twisted and the graph embeds into 
a copy of the projective sphere. 

We now proceed with the induction step. One notes that by independence of 
the pairs at different edges, the left hand side of (|3.2p factors into the product 
corresponding to connected components of F. It is therefore enough to consider 
connected F. 

If F has two vertices that are joined by an edge, we can use cyclicity of to 
move the variables that label the edge to the first positions in their cycles, say Xi 
and Xj^^i and use (|2.6p or (|2.7p to eliminate this pair from the product. The use 
of relation (|2.6p is just that of gluing the two vertices together removing the edge 
X which is labeled by the two appearances of Z. Relation (|2.7|) glues together the 
two vertices, removing the edge x, and the reversal of orientation across the edge 
is given by the conjugate (see Figure [2]). These geometric operations reduce n and 
m by one without changing the Euler characteristic: the number of edges and the 
number of vertices arc reduced by 1; the faces are preserved - in the case of edge 
flip in Fig. [21 the edges of the face from which we remove the edge, after reduction 
follow the same order. 

Therefore we will only need to prove the result for the single vertex case of the 
induction step. 

We wish to show that 

(3.3) E{X,X2X3X4 . . . X2n) = (-2)'^(i^)4"-i , 

where F is a one vertex Mobius graph with arrows (half edges) labeled by Xk- We 
will do this by induction, there are two cases: 
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Figure 2. A Mobius graph with two vertices connected by a rib- 
bon may be reduced to a Mobius graph with one less vertex and 
one less edge. In these graphs the ". . ." are to mean that there 
are arbitrary numbers of other edges at the vertex, and the edges 
drawn are connected to other vertices. The top graph is an exam- 
ple of this reduction when the connecting ribbon is untwisted, in 
this case the two vertices are glued together with no other changes 
in the ribbons. The bottom graph is an example of this reduction 
when the connecting ribbon is twisted, in this case the two ver- 
tices are glued together and the order and orientation (twisted or 
untwisted) of the ribbons on one side are reversed. 



Case 1: Xi= Xj for 1< j < 2n, 

E{XiX2 . . . + i . . . X2n) = E(XiXi)E(5R(X2 . . . + i . . . X2n)) 

(3.4) = 4E(3ff(X2 . . . X,_i)5ft(X,+i . . . X2n)) ■ 

This corresponds to the reduction of the Mobius graph pictured in Figure 
[31 which splits the single vertex into two vertices. The Eulcr characteristic 
becomes X2 = 2 — (n — 1) + /i = xi + 2. By the induction assumption we 
find 



E(Xi . . . X2n) = 4 



(-2)^2 = 4""2(-2)'^i+2 = 4""i(-2)' 
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A 

•* 

Figure 3. Here we have an untwisted ribbon of the Mobius graph 
returning to the same vertex, this edge is removed in our reduction 
procedure giving us a Mobius graph with one more vertex and one 
less edge. 

A 

Figure 4. Here we have a twisted ribbon of the Mobius graph 
returning to the same vertex, this edge is removed in our reduction 
procedure giving us a Mobius graph with one less edge and giving 
a reverse in both the order and orientation (twisted or untwisted) 
of the ribbons on one side of the removed ribbon. 





Case 2: Xi = Xj for 1 < j < 2n, 

E(XiX2 . . . + i . . . X2n) = E(XiXi)E(X,_i . . . M. + i . . . X2n) 

(3.5) = (-2)E(X,_i . . . Wj + i . . . X2n) ■ 

This corresponds to the reduction of the Mobius graph pictured in Figure 
Ul which keeps the single vertex and flips the order and orientation of the 
edges between Xi and Xj. The Eulcr characteristic becomes X2 = 1 — (n — 
1) + /i = Xi + 1- By the induction assumption wc find 

E(Xi . . . X2n) = (-2) W''-^'^-^{-2YA = (-2)-i4"-i(-2)>^i+i = 4"-i(-2)xi . 

Note that taking away an oriented ribbon creates a new vertex. The remaining 
graph might still be connected, or it may split into two components. If taking 
away a loop makes the graph disconnected, then the counts of changes to edges and 
vertices are still the same. But the faces need to be counted as follows: the inner 
face of the removed edge becomes the outside face of one component, and the outer 
face at the removed edge becomes the outer face of the other component. Thus the 
counting of faces is not affected by whether the graph is connected. 
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With these two cases checked, by the induction hypothesis, the proof is com- 
pleted. 

□ 

3.3. Bipartite Mobius graphs and quaternionic Gaussian moments. To 

deal with quaternionic Wishart random matrices, we need to consider a special 
subclass of quaternionic Gaussian variables from Theorem 13.11 Suppose the 2n- 
tuplc {X±i, X±2, . . . , X±n) consists of n pairs of random variables taken with repe- 
tition, from the set {Zi, Z2, . . . , Zn} of independent quaternionic Gaussian random 
variables. Note that in contrast to the setup for Theorem 13.11 here all Z's are 
without conjugation. Fix a sequence ji, J2, • ■ • ,jm of natural numbers such that 
jl + ■ ■ ■ + jm = n. Theorem 13.11 then says that 

E(JR.{X_iXi ■ ■ ■X_jj^Xj^) X 3fJ(X_jj_iXjj_|_i • • ■ X^j^^j^Xj^+j^) X ■ • ■ 

• ■ • X • • ■ X^nXn)) 

where T is the Mobius graph with m vertices and with edges labeled by 
X±i, . . . ,X±n that describes the pairings between the variables under the expec- 
tation. Our goal is to show that the same formula holds true for another graph, a 
bipartite Mobius graph whose edges are labeled by n pairs {X-j,Xj), 1 < j < n. 

The bipartite Mobius graph has two types of vertices: black vertices and white 
(or later, colored) vertices with ribbons that can only connect a black vertex to a 
white vertex. (As previously, the ribbons may carry a "flip" of orientation which 
we represent graphically as a twist.) To define this graph, we need to introduce 
three pair partitions on the set {±1, . . . , ±n}. The first partition, 6, pairs j with 
—j. The second partition, cr, describes the placement of 5R: its pairs are 

{l,-2},{2,-3},...,{ji,-l}, 

{ji + 1, -i:h + 2)}, . . . , {jl + J2, -Ui + 1)}, 



rn — 1 m — 1 

Jfe- -(2 + E ^■''Ol. . . . , K -(1 + E Jfe)}- 

fc=l k=l k=l 

The third partition, 7, describes the choices of pairs from Zi, . . . , Z„. Thus {j, k} G 
7 if Xj = Xk when jk > or Xj = Xk if jk < 0. 

We will also represent these pair partitions as graphs with vertices arranged in 
two rows, and with the edges drawn between the vertices in each pair of a partition. 
Thus 

12 n 
mm m 

(3.6) I I I 

• • • 

-1 -2 ■ ■ ■ -n 

and 

1 2 3 ... n Jl+l 31+2 ... ji+32 jl+32 + 1 
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Figure 5 . Representation of a bipartite Mobius graph, the edges 
drawn arc ribbons that are either twisted or untwisted. 



Consider the 2-regular graphs S U j and 5 U a. We orient the cycles of these 
graphs by ordering on the left-most vertical edge of the cycle. For example, 




We now define the bipartite Mobius graph by assigning black vertices to the m 
cycles of (5 U cr, and white vertices to the cycles of (5 U 7. 

Each black vertex is oriented counter-clockwise. For each black vertex we follow 
the cycle of 5 U cr, drawing a labeled line for each clement of the partition. The 
lines corresponding to are adjacent, and will eventually become two edges of 
a ribbon. 

Each white vertex is oriented, say, clockwise; we identify the graphs that differ 
only by a choice of orientation at some of the white vertices. For each white vertex, 
we follow the corresponding cycle of (5 U 7, drawing a labeled line for each clement 
of the partition. The lines corresponding to are adjacent, but may appear in 
two different orders depending on the orientation of the corresponding edge of S on 
the cycle. 

The final step is to connect pairs (j, —j) on the black vertices with the same pairs 
on the white vertices. This creates the ribbons, which carry a flip if the orientation 
of the two lines on the black vertex docs not match the orientation of the same 
edges on the white vertex. 

The individual edges pictured in Figure[5]are ribbons and are labeled as in Figure 
[6l We allow twists of ribbons to propagate through a white vertex calling the two 
bipartite Mobius graphs equivalent. 

Suppose the 2n-tuple {X±i, X±2, . . . ,X±n) consists of random variables taken 
with repetition, from the set {Z'l, Z2, . . . , Z„} of independent quaternionic Gaussian 
random variables. Let A4 = X±2, ■ • ■ , X±n) denote the set of all bipartite 

Mobius graphs F that correspond to various ways of pairing all repeated Z's in the 
sequence {X±i, X±2, . . . ,^±n); the pairs are given by adjacent half edges at each 
white vertex. (See the preceding construction.) = if there is a Zj that is 
repeated an odd number of times. 
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X 



X. 



X. 



x_ 



X 



X, 



x_ 



Figure 6. Example of the labeling we use for the edges emanating 
from a black vertex in a bipartite Mobius graph. 



Theorem 3.2. 

• ■ • X K(X_j^_j2 j„i-i-l ■ ■ ■ X^nXn)) 



(The right hand side is interpreted as when A4 ~ 9.) 



Proof. The proof is fundamentally the same as that of the Wigner version of this 
theorem. In view of (|2.9p and (|2.10p . it suffices to consider {X±i, . . .X±n} that 
form n independent pairs, and show that 

(3.7) E(5R(X_iXi • • • X^,,X,J X 5R(X_,,_iX,,+i • • • x ■ ■ ■ 

where F is the bipartite Mobius graph that describes the pairings. 

We will prove (|3.7p by induction; to that end we first check that with n = 1, m = 
1 we have E{XX) = (-2)^ in agreement with ([XT)) . 

If F has two black vertices connected together by edges adjacent at a white 
vertex, we can use the cyclicity of to move the variables that label the respective 
edges and share the same face to the first position in their cycles, so that we may 
call them X_i and either Xj or X^j. We now use relations (|2.6p and (j2.7p to 
eliminate the pair from the product: 

E (5R(X_i • • • X,,MX, . . . X,,+,,X.,)) = E(Xi • • • X,,+,,X.,) , 



E (5R(X_i . . . X,,)n{X, . . . = E(X,,X_,, . . . X,X, ■ ■ ■ . 

The use of relation (|2.6p corresponds to that of gluing together the two ribbons along 
the halves adjacent at the white vertex, and gluing together the corresponding black 
vertices (see Figure [7]). The use of relation (|2.7p corresponds to the same gluing, 
but in this case one of the ribbons has an orientation reversal in it, resulting in 
an orientation reversal for the remaining sides (see Figure [8]). These geometric 
operations reduce n and m by one without changing the Euler characteristic: both 
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Figure 7. Here two black vertices arc connected together through 
untwisted edges adjacent at a white vertex. This bipartite Mobius 
graph reduces to one with one less vertex and one less edge. The 
reduction is found by gluing the two black vertices together and 
gluing the two ribbons together along their adjacent sides, here 
labeled by X-i and Xj. The same reduction would apply if the 
two edges were twisted as we could pass this twist through the 
white vertex. 




Figure 8. Here two black vertices are connected together through 
one twisted edge and one untwisted edge adjacent at a white vertex. 
This bipartite Mobius graph reduces to one with one less vertex 
and one less edge. The reduction is found by gluing the two black 
vertices together, gluing the two ribbons together along the sides 
adjacent at the white vertex, and reversing both the order and 
orientations of remaining ribbons on the second black vertex. 



the number of edges and the number of vertices are reduced by one while the number 
of faces is preserved. 

Therefore we will only need to prove the result for the single black vertex case 
of the induction step. We wish to show that 

(3.8) E(X_iXiX_2^2---^-„^„) = 4"-i(-2)'^(i^) , 

where F is a bipartite Mobius graph with a single black vertex and half ribbons 
labeled by X±k. We will do this by induction, there are two cases: 
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Figure 9. Here we have a black vertex with two ribbons, both 
twisted or both untwisted, adjacent at the same white vertex. The 
reduction glues these two ribbons together along their common 
side. The result is a bipartite Mobius graph with one more ver- 
tex and one less edge. The resulting graph may or may not be 
disconnected at this point. 



Case 1: X_i = Xj for 1 < j < n, 

E(X_iXi • • • X^,X, • • ■ X„) = E(X_iX,)E(K(Xi • ■ • ■ • • X„)) 

= 4E(3?(X_,Xi ■ • •X,_i)5R(X_,_iX,+i • • •X„)) . 

This corresponds to the reduction of the bipartite Mobius graph pictured 
in Figure [9l which for j > 1 splits the single black vertex into two black 
vertices, and glues the two edges labeled as X-i and Xj together. 

The edges (1,— — j) are adjacent at the white vertex and appear 
either in this order, or in the reverse order. Thus after the removal of 
{ — 1, j} we get an ordered pair of labels (1, — j) or (— j, 1) to glue back into 
a ribbon. On the black vertex, due to our conventions the edges of the 
ribbons appear in the following order 

((-1, 1), (-2, 2), . . . , j), (-fci, fci), . . . , i-kr, kr)). 

Once we split the black vertex into two vertices with the edges of rib- 
bons given by ((-1, 1), (-2, 2), . . . , (-j, j)) and ((-fci, fci), . . . , (-fc^, fc^)), 
the removal of { — creates a new pair (1,— j) which wc use to create 
the ribbon to the white vertex. [After this step, we relabel all edges to use 
again ±1, ±2, . . . consecutively.] 

We note that the number of faces of the new graph is the same as the 
previous graph - the face with the sequence of edges 

(■ • • , Xk^,X^i,Xj,X^ki , • • • ) 

becomes the face with edges (. . . , Xk^ , X^ki , ■ • ■ ) on the new graph. The 
Euler characteristic becomes X2 = {vi + 1) — (ei — 1) -|- /i = Xi + 2 where 
vi,ei and /i are the number of vertices, edges, and faces of F. By the 
induction assumption we then find 

E(X_iXi---X„) =4 f4("-i)-2(_2)X2l ^4«-2(_2)xi+2 =4"-i(_2)xi . 
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Figure 10. Here we have a black vertex with two ribbons, one 
twisted and the other untwisted, adjacent at the same white vertex. 
The reduction glues these two ribbons together along their common 
side. The result is a bipartite Mobius graph with one less edge, 
and with the ribbons on one side of the removed ribbon now with 
reversed order and orientations. 



Case 2: X_i = X^j for 1 < j < n, 

EiX^iXi---X.,X,---Xn) 

= E(X_iX_j)E(Xj_iX_j+i • • • XiXjX-j-iXj+i ■ ■ ■ Xn) 

= (— 2)E(Xj^iX_j+i • • • XiXjX-j^iXjj^i ■ ■ ■ Xn) ■ 

This corresponds to the reduction of the Mobius graph pictured in Figure 
I10| which switches the order and the orientations of the ribbons on one 
side of the black vertex from the ±1 and ±j ribbons and glues the two 
edges adjacent at the white vertex together as shown. As previously, the 
removed edges are adjacent at the white vertex. At the black vertex, the 
labeled lines for the construction of the bipartite graph change from the 
sequence 

(-1, 1), (-2, 2), . . . , (-J + 1, J - 1), (-J, j), (-fci, fci), . . . , {~kr, kr) 
to the sequence 

(j - 1, -J + 1), . . . , (2, -2), (1, j), (-fci, fci), . . . , {-K, kr) 

which then needs to be relabeled to use ±1, . . . , in. Again the number of 
faces on the bipartite graphs is preserved: the face with edges 

(...,2,-1, -J, /ci,...) 

becomes the face 

(...,2,1, J, fci,...). 

The Euler characteristic becomes X2 = t'l ~ (ei — 1) + /i = Xi + 1- By the 
induction assumption we then find 

E(X_iXi ■■■Xn) = (-2) [4("-i)-i(-2)x=l = (-2)4"-2(-2)xi+i = 4"-i(-2)xi . 



With these two cases checked, by the induction hypothesis, the proof is com- 
pleted. 

□ 
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One should note that this is fundamentally the same proof as in the Wigner 
case, however in this case the geometric reduction is given by gluing together two 
ribbons, while in the Wigner case the geometric reduction is the elimination of one 
ribbon at a time. The inductive steps remain the same. 



4. Duality between real and symplectic ensembles 

By A4MxAf(IHI) we denote the set of all AI x N matrices with entries from H. For 
A e MmxnC^), the adjoint matrix is A*j := Aj^i. The trace is tr(A) = J2f=i ^jj- 
Since the traces tr(A) may fail to commute, in the formulas we will use 5R(tr(A)), 
compare to [TU] • 

4.1. Duality between GOE and GSE ensembles. The Gaussian orthogonal 
ensembles consist of square symmetric matrices, Z, whose independent entries are 
independent (real) Gaussian random variables; the off diagonal entries have variance 
1 while the diagonal entries have variance 2. One may show that 

Theorem B. For Z from the N x N Gaussian orthogonal ensemble: 

^^E(tr(Z^-^)tr(Z-'-^) . . .tr(Z-'")) = ^ TV^^^^) , 

r 

where the sum is over labeled Mobius graphs T with m vertices of degree 
Jij J2, • ■ • , Jm, x(r) is the Euler characteristic and ji+j2 + - ■ ■+jm = 2n. More gen- 
erally, i/Zi, . . . , Zs are independent N x N GOE ensembles and t : {1, 2, . . . , 2n} 
{1, . . . , s} is fixed, then 



'T:(tr(Zt(i) . . . Zt(/3^))tr(Zt(„2) • ■ • Zt(/32)) x 



(r) 



•••xtr(Z,(„„)...Z,(^„)))=^iV>^(r 

r 

where ai = 1 , aj,. = ji + J 2 H h jfe-i + 1, and = ji + ^2 H \~ jk denote the 

ranges under the traces, and where the sum is over labeled color-preserving Mobius 
graphs T with vertices of degree ji, j2, ■ ■ ■ ,jm whose edges are colored by the mapping 
t. If there are no T that are consistent with the coloring we interpret the sum as 
being 0. 

The single color version of this Theorem was given in Ref. [5] . 

The Gaussian symplectic ensembles (GSE) consist of square self-adjoint matrices 

(4.1) Z = [Z,,,] , 

where {Zi^ : z < j} is a family of independent H-Gaussian random variables (|2.3p . 
{Zi i} is a family of independent real Gaussian random variables of variance 2, and 
Zij = Zj^i. The law of such a matrix has a density Cexp(— i5Rtr(x^)), which is 
supported on the self-adjoint subset of Mnxn{^)- {C — C{N) is a normalizing 
constant.) 

In this section we will demonstrate a multi-matrix version of the theorem of 
Mulase and Waldron [21] from our Wick formula. This represents an improvement 
to the argument of Mulase and Waldron as we will not need to rely on labelings of 
the vertices by the quaternions 1, i, j, and k. This part of their argument is encoded 
in relations (|2.4ll2.7p . 
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We will compute here the expected values of the real traces of powers of a 
quaternionic self-dual matrix in the Gaussian symplectic ensemble (|4.ip where the 
off-diagonal matrix entries are (|2.3p . 

The basic theorem we will prove is 

Theorem 4.1. For Z from the N x N Gaussian symplectic ensemble: 

where the sum is over labeled Mobius graphs T with m vertices of degree 
ji, j2, . . . ,jm, is the Euler characteristic and ji + .72 + ' ' ' + jm = 2n. 

More generally, suppose Zi, . . . , are independent N x N GSE ensembles and 

t : {I,..., 2n} ^ {1, . . . , s} is fixed. Let ai = 1, at = ji + j-i ^ h jk-i + 1, 

Pk = ji + j2 + • ■ ■ + jk denote the ranges under traces. Then 

(4-3) (4^|„-.„ iE(^(tKZt(i) • . . Z,(^,)))5?(tr(Z,(„,) . . . Z,(^,))) x . . . 

• • • X 5R(tr(Z,(„„) . . . Zt(p^)))) = ^(-27V)>^(r) ^ 

r 

where the sum is over labeled color-preserving Mobius graphs V with vertices of 
degree ji, j2, . . . , jm that are colored with s colors by the mapping t. As previously, 
x(r) is the Euler characteristic and ji + j2 + ■ ■ ■ + jm ~ 2n. (If there are no T that 
are consistent with the coloring t, we interpret the sum as Q.) 

In view of Theorem[Bl Theorem 14.11 gives the duality between the GOE and GSE 
ensembles of random matrices as ^ —2N. 

Example 4.1. To illustrate the multi-matrix aspect of the theorem, suppose Zi, Z2 
are independent N x N GSE ensembles and Xi , X2 are independent N x N GOE 
ensembles. For p,r,k > let bp,r,k{N) := E((tr(Xi^Xf ))'=). By Theorem [H 

bp,r,k{^) = N'^P+''^^^''J2r^^^^^- Therefore Theorem O implies that the mo- 
ments for the independent GSE ensembles are determined from the corresponding 
moments of independent GOE ensembles by the "dual formula" 

E((tr(Z?Pzr))''0 = (-2)(P+'-i"^6p,,,,(-27V). 

In particular, 

E(tr(Z?^)) = {~2r~\,o.AN). 

This is equivalent to [HI Theorem 6] with Cp{N) = E(tr(Z2p))/2P; an extra factor 
of 2^ in our formula is due to the fact that the symplectic ensemble in [T7] has the 
variance of 1/2 instead of our choice of 1. 
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Proof of Theorem \4-l\ We begin by expanding out the traces in terms of the matrix 
entries 



(4.4) 

I 

i 



^E(5R(tr(Z^-^))5R(tr(Z^-^)) . . . 5ft(tr(Z^"))) 



(4n)"-" ^(^(^ai.a2^a2,a3 ■ • ■ Zaj^ ,ai)^{Zbi ,b2 ■ ■ ■ Zbj^,bi) X . . 

^ •••x5R(Z,„,,...Z,^,„,,J), 
1 < fli, a2, . . . , < Jy 

1 < bi,b2,...,b,, < N 



1 < Ci, C2, . . . , Cj^ < N 
Note that essentially the same expansion applies to (|4.3p . except that the consec- 



utive entries ^ must now be labeled also by the "color" t. It is then better to 



index the products by the cycles of a permutation. Put 

cr = (1, . . . ,/3i)(a2, . . . ,/32) • • • (ttm, . . . ,2ri,) . 
In this notation, the expansion of the left hand side of (|4.3p takes the form 

E E n^(n[z*o)]ao),aM.))) 

o:{1...2n}-^{l...Ar} \cei7 jec 

(Since has the cyclic propertv l2.21 this expression is well defined.) 

From (|2.10p it follows that the right hand side of (|4.4p can be expanded as a 
sum over all pairings, and we can assume that the pairs are independent. 

Of course, pairings that match two independent random variables do not con- 
tribute to the sum. The pairings that contribute to the sum are of three different 
types: pairs that match Z with another Z at a different position in the product, 
pairs that match Z with Z, and pairings that match the diagonal (real) entries. 

We first dispense with the pairings that match a diagonal entry Zi^i = ^ with 
another diagonal entry Zj,j = ^. Since real numbers commute with quaternions, 

^{EiqoCqi^q2)) = EieMqoqiq.J = 25R(gogi92). 

On the other hand, using the cyclic property of 5R(-) and adding formulas (|2.4p and 
(|2.5p . we see that the same answer arises from 

E^{qoZqiZq2)+m{qoZqiZq2) = E^{Z qi Z q2qo) + m{Z qi Z q2qo) = 23fJ(gig29o). 

So the contribution of each such diagonal pairing is the same as that of two matches 
of quaternionic entries {Z,Z) and {Z,Z). 

Thus, once we replace all the real entries that came from the diagonal entries 
by the corresponding quaternion-Gaussian pairs, we get the sum over all possible 
pairs of matches of the first two types only. We label all such pairings by Mobius 
graphs, with the interpretation that pairings or random variables {Z, Z) correspond 
to twisted ribbons, while pairings {Z, Z) correspond to ribbons without a twist. 
The pairs of variables at different ribbons can now be assumed independent. In the 
multi-matrix case, the edges of the graph at each vertex are colored according to 
the function t, which restricts the number of available pairings. 
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Figure 11. For the Wishart applications we will label the sides 
of the ribbons emanating from a black vertex. One side receives 
the label Xjn,mk while the other receives the label Xm^uk+i- Note 
that the interior index, mk, of the two labels on each ribbon match; 
while the exterior index, Uk, of the labels on ribbons adjacent at a 
black vertex match. 



We now relabel the by Zak,ak+i = Xk, ^h^.b^+i = Xj^+k, -^€^,0^+1 = 
Xjj^+j24_...j^_j+fc. Our claim is then that given a Mobius graph F with vertices of 
degree ji , j2 , ■ • ■ , and jm with edges labeled by X^ , satisfies 

(4-5) (4^|„_„, E(3^(Xi . . . . . . . ..X2n)) 

= (-27V)^(r)7v-/(r) ^ 

where /(F) is the number of faces of F. The N~f^^'> terms are removed by the 
summations in ()4.4|) . the relations given by the edges of F reduce the number of 
summations leaving us with /(F) sums from 1 to A^. Collecting the powers of N 
we find that this is the same as p.2p . □ 

4.2. Duality between the Laguerre Ensembles. Wc will now consider the anal- 
ogous problem for quatcrnionic Wishart random matrices. These are N x N ma- 
trices of the form 

(4.6) W = X*X, 

where X is an M x N matrix with independent H-Gaussian entries (|2.3p . The Wick 
formula (|2.10[) in this case is expressed in terms of bipartite Mobius graphs defined 
in section [3. 31 although in this case it is helpful to view our labeling of the ribbons 
around each black vertex in a slightly different way (see Figure [TT|) . Our heuristic is 
the following: when edges are glued to a white vertex we identify the labels on the 
interior of each ribbon as all equal, while the labels on the outside of each ribbon 
are equal to the others which border the same face. The reason for this dichotomy 
is that X is not a square matrix and therefore the rij's and rrij's arc to be treated 
as separate families of variables. 

In the multi-matrix case, we color the ribbons by the number of a matrix, and 
color the white vertices so that only the ribbons of matching color can be attached. 
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X 



niiTii 



X = Hi 



O 



Figure 12. The only possible bipartite Mobius graph with a sin- 
gle black vertex of degree one. In this case the graph is embedded 
on a copy of the Riemann sphere. 



M N 



For example, 

E(tr(W)) = E^(^n™^"-) = 

m— 1 n—1 

corresponds to the bipartite Mobius graph in Figure [T21 As another example 

M N 

E(tr(w2))= j2 E E(^:,„Hx.un2^:.™.^™.nj 

mi, 7712 — 1 711,712 — 1 

M N 

= E E [l'^("'^ni77ii-'^mi7i2)I^(-'^n2m2^™2ni) 

7711, 7712 — 1 771,772 — 1 

+ E(X„^ ^77 2 7712 )-'^(^niin2 ^7712)11 ) 

= 16NM^ - 8NM + 16N^M , 

corresponds to the bipartite Mobius graphs in Figure [T31 
In this case our theorem takes the form. 

Theorem 4.2. With M = XN we find 

E(tr(WJi)tr(WJ'^) • • •tr(W^'")) = 7V"-"4"-'" ^ A"'^^) (-2iV)'^(^) 



where the sum is over the bipartite Mobius graphs with black vertices having degrees 
jii j2, ■ ■ • , jm, n = ji + j2 + ■ ■ ■ + jm *s the number of edges, w(T) is the number of 
white vertices, and x(r) is the Euler characteristic. 

More generally, suppose that Wi,W2, . . . , Ws are N x N independent quater- 
nionic Wishart with parameters Mi,...,Ms. Denote Xj = Mj/N. Fix t : 
{!,..., n} -> {!,..., s}, and let ai = l,ak = ji + ^2 + • • • + jfc-i + 1, 
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ni = n2 mi^rrij 




fifiinj ni = ni n2 = n2 mi = m2 



Figure 13. The three possible bipartite Mobius graphs with a 
single black vertex of degree two. In the top two cases the graph is 
embedded on a copy of the Ricmann sphere, while in the bottom 
case the graph is embedded on a copy of the projective sphere. 



/3fe = ji + h 
1 



ik denote the ranges under traces. Then 



(4-7) (4^)„^.„ 'E(^(tKW,(i) . . . W,(^,)))jR(tr(W,(,,) . . . W,(fe))) x 



X 5R(tr(W,(„,„) . . . W,(^„^)))) = ^ Ar^(^^A-^(^) 



.A 



{~2N) 



xiT) 



where the sum is over the bipartite Mobius graphs with m black vertices of degrees 
jii ■ ■ ■ T jm, whose ribbons are colored according to t by one of the colors 1, . . . , and 
with white vertices that are colored to match the ribbons. Here Wj{T) is the number 
of white vertices of color j , x(r) is the Euler characteristic and ji+j2 + - ■ ■+.jm = n 
is the number of edges /ribbons. 

Proof. We begin by expanding out the traces in terms of the matrix entries of X 
where W = X*X and X is an M x N matrix of quaternionic Gaussian random 
variables 
(4.8) 

Arm — n 

E(tr(WJi)tr(W^"^) • • • tr(W^")) = 

^ ^„_,„ E(jR([X*]a^Ai [XJ^iaa [X*]a2A2 WAaaa ■ ■ ■ [^*]aj^ Aj^[X]a 



1 < ai, 02, . . . , aj^ < 
1 < Ai,A2,...,Aj, < M 

1 < 61,62,..., <iV 
1 < Bi,B2,...,Bj, < M 



5ft([X*]b,B,...[X]B,,6j...K([X*],,c,...[X: 



c,-, 



1 < Ci,C2,...,Q„ < M 

Note that a similar expansion would be found for a multi matrix expression, see 
below. 

From (|2?T0l) it follows that (|48l) can be expanded as a sum over all pairings, and 
we can assume that the pairs are independent. 
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We label all such pairings by bipartite Mobius graphs, that is Mobius graphs 
with labeled black vertices of degree ji , j2 , . . . , jm , and unlabeled white vertices of 
any degree. The half edges of each graph arc labeled by XA^ai = ^-i, -^Aia2 = 
and so on. The relations given by this Mobius graph reduce the number of sums 
over capital indices (which go from 1 to M) to wiT) the number of white vertices, 
and reduce the number of sums over lower case indices (which go from 1 to N) to 
/(r) the number of faces of F. Therefore from Theorem 13.21 with M = XN, the 
total contribution from the Wishart Mobius graph F to the sum is 

as x(r) = m + w{T) -n + f{T). 

The multimatrix version of this proof requires no conceptual changes, but the 
notation becomes more cumbersome. In this case it is more convenient to index 
the products by the cycles of a permutation. Put 

(T = (1, . . . ,/3i)(a2, . . . ,/32) • . • {am, ■■■,n). 

In this notation, the expansion on the right hand side of (|4.8p is replaced by 



E E n^(n[Wta)]ao),a(.o-))) 

a:{l...ri}^{l...Af} \cecr j<£c 

E E n^'^n[^*(i)]60)a(j)[^tO)]Kj),^>('^0))) 

a:{l. .. n}^{l. .. N} beB{t) \ce<J jGc 

where 6(t) -{!,..., Mt(i)} x {1, . . . , Mt^^)} x • ■ • x {1, . . . , Mt(„)}. 

From (|2.10p the expected value can again be expanded as a sum over the pair- 
ings, where we can assume that the pairs are independent random variables and 
the pairings connect only the random variables of the same color t{j). For each 
bipartite Mobius graph F, the sum over b E B{t) contributes a factor of MJ''"^^ per 
color, where Wj{T) is the number of vertices of color j in F, while the sum over a 
contributes N-f^^\ as in the one-matrix case. Therefore from Theorem 13.21 the left 
hand side of (14.71) becomes 



r 

s 

= 7Vm-n Y Ar«'(r)+/(r) (_2)x(r) "Q xj^ 



= E 7v"'(^)+™+-''(^^-"(-2)'^(^) 11-^7'^^^ 

S S 

r j=i r 3=1 



□ 
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We remark that the result is again in duahty to the formula for real Wishart 
matrices, where the analogue of (|4.7|) has the following form: 

(4-9) ^^lE(tr(W,(i) . . . W,(^,))tr(W4(„,) . . . Wi(^,)) x . . . 

... X tr(W,(„„) . . . W,(,,„))) = ^ Ar^^^Ar^'^ . . . A-(i^)7V>^(r). 

r 

For the one-matrix case this can be read out from |10| . For the multivariate case 
this is a reinterpretation of [H Theorem 2.10]. 

The following illustrates a one-matrix version of this duality. 

Example 4.2. Suppose Wjj is x iV real Wishart matrix with parameter M, and 
let 

Ql_^jM,N) = E(tr(W;^)tr(W^^) . . .tr(W^'")). 

If Wh is N xN quaternion Wishart matrix with parameter M, and n = + - ■ ■ +jm , 
then ((^ and (|i?7)) give 

E(tr(W£)tr(W^=)...tr(W^")) = (-2r-'"Q« (-2M, -2iV). 
This is equivalent to [10[ Corollary 4.2] , as with X ~ [ji, j2, ■ ■ ■ , jm) in their notation 

Ql_,jM, N) = 2-"E(tr(W^^)tr(W^^) . . . tr(W>r )) 

= (-l)"-™2-™Q«^...^^-^^(-2M,-27V); 

an extra factor of 2" in our formula is due to the fact that the quaternion Gaussian 
law in [TOj has the variance of 1/2 instead of our choice of 1. 

Acknowledgements. The authors thank Gerard Letac for citation [T^]. We would 
like to thank the referee for several substantial corrections. 
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